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Abstract 

We present a unified global analysis of neutrino oscillation data within the framework of the 
four-neutrino mass schemes (3-|-l) and (2-|-2). We include all data from solar and atmospheric 
neutrino experiments, as well as information from short-baseline experiments including LSND. If 
we combine only solar and atmospheric neutrino data, (3-1-1) schemes are clearly preferred, whereas 
short-baseline data in combination with atmospheric data prefers (2-1-2) models. When combining 
all data in a global analysis the (3-1-1) mass scheme gives a slightly better fit than the (2-1-2) 
case, though all four-neutrino schemes are presently acceptable. The LSND result disfavors the 
three-active neutrino scenario with only ArngQ^ at 99.9% C.L. with respect to the four- 

neutrino best fit model. We perform a detailed analysis of the goodness of fit to identify which 
sub-set of the data is in disagreement with the best fit solution in a given mass scheme. 
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I. INTRODUCTION 


From the long-standing solar and atmospheric BBI-B neutrino anomalies 

we now have compelling evidence that an extension of the Standard Model of particle physics 
is necessary in the lepton sector. The most natural explanation of these experiments is 
provided by neutrino oscillations induced by neutrino masses and mixing with neutrino mass- 
squared differences of the order of eV^ and ~ 3 x 10“^ eV^. Explaining 

also the evidence of > Vg oscillations with a mass-squared difference AtUlsno ~ 1 
reported by the LSND experiment |]^ requires an even more radical modihcation of the 
Standard Model. Currently this experiment is left out in most analyses of neutrino data. At 
the moment the LSND result is not conhrmed nor ruled out by any other experiment, and 
therefore it is reasonable to see more quantitatively its impact on the physics of the lepton 
sector. 

If all the three anomalies are explained by neutrino oscillations, and the possibility of CPT 
violation is neglected |^|, we need at least four neutrinos to obtain the three required mass- 
squared differences. In view of the LEP results the fourth neutrino must not couple to the 
Z-boson. Such a sterile neutrino with a mass in the electron-volt range has been postulated 
originally to provide some hot dark matter suggested by early CODE results |^, |^, p!3|. 


and after the LSND result many four-neutrino models have been proposed [|T^, |T9 

A quite complete list on four-neutrino references can be found at . 

A very important issue in the context of four-neutrino scenarios is the question of the 
four-neutrino mass spectrum. Two very different classes of four-neutrino mass spectra can 
be identihed. The hrst class contains four types and consists of spectra where three neutrino 
masses are clustered together, whereas the fourth mass is separated from the cluster by the 
mass gap needed to reproduce the LSND result. The second class has two types where 
one pair of nearly degenerate masses is separated by the LSND gap from the two lightest 
neutrinos. These two classes are referred to as (3-1-1) and (2-1-2) neutrino mass spectra, 
respectively |^. All possible four-neutrino mass spectra are shown in Fig. |l[ 

One important theoretical issue in these models is how to account for the lightness of 
the sterile neutrino which, ordinarily, should have mass well above the weak scale. The 
simplest possibility is to appeal to an underlying protecting symmetry, getting, moreover, 
the LSND mass at one-loop order only |]T^ Alternatively, the lightness of the sterile 


neutrino may follow from volume suppression in models based on extra dimensions jlTl , [18 


As for the maximal atmospheric mixing angle, it follows naturally in the models of Refs. O, 


1^, h8| since to hrst approximation the heaviest neutrinos form a quasi-Dirac pair whose 


components mix maximally. Finally the splittings which generate solar and atmospheric 
oscillations arise due to breaking of the original symmetry (for example due to additional 
loop suppression) |^, or due to R-parity breaking These models lead to a (2-1-2) 
scheme. 

One important feature of (3-1-1) mass spectra is that they include the three-active neu- 
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trino scenario as limiting case. In this case solar and atmospheric neutrino oscillations are 
explained by active neutrino oscillations, with mass-squared differences Am^QL and 
and the fourth neutrino state gets completely decoupled. We will refer to this scenario as 
(3-1-0). The (3-1-1) scheme can be considered as a perturbation of the (3-1-0) case: a small 
mixture of and with the separated mass state can account for the oscillations observed 
by LSND. In contrast, the (2-1-2) spectrum is intrinsically different from the three-active 
neutrino case. A very important prediction of this mass spectrum is that there has to be 
a significant contribution of the sterile neutrino either in solar or in atmospheric neutrino 
oscillations or in both. More precisely, in the (2-1-2) case the fractions of sterile neutrino 
participating in solar and in atmospheric oscillations have to add up to one . 


Based on semi-quantitative arguments it has been realized for some time [23, M, 25, 26 


that it is difficult to explain the LSND result in the framework of (3-1-1) schemes because 
of strong bounds from negative neutrino oscillation searches in short-baseline (SBL) exper¬ 
iments, and therefore the (2-1-2) scheme was considered as the preferred one. Recent exper¬ 
imental developments lead to a renaissance of the (3-1-1) mass schemes ^ |^. First, 
a new LSND analysis (see last reference of [^) resulted in a shift of the region allowed 
by LSND to slightly smaller values of Am^sNO) which makes the (3-1-1) schemes somewhat 


less disfavored. However, in Refs. |]28| , ^ it was shown within a well defined statistical 
analysis that a bound implied by SBL experiments is in disagreement even with the new 
LSND allowed region at the 95% C.L. in (3-1-1) schemes. Second, the high statistics data 
from Super-Kamiokande started to exclude two-neutrino oscillations into a sterile neutrino 
for both solar as well as atmospheric neutrinos |^, which constitutes a problem for (2-1-2) 
mass schemes. Concerning the solar data, the trend to disfavor oscillations into a sterile 
neutrino recently became supported by the beautiful result of the SNO experiment P, 
However, a unified analysis of solar and atmospheric neutrino data performed in Refs. |^, |3^ 
showed that the goodness of fit of the (2-1-2) mass scheme is still acceptable. 

In this work we perform for the first time a global analysis of all the relevant neutrino 

oscillation data in a four-neutrino framework. We will use the fit of the global solar neutrino 

data presented in Ref. |^, which includes Super-Kamiokande ||l[|, Homestake @], SAGE [Q, 

GALLEX and GNO [Q and SNO P|. Further we include data from the atmospheric neu- 

(-) 


trino experiments Super-Kamiokande [0 and MACRO |p, data from the SBL u] 


(-) 
l^e 

appearance experiments LSND [^, KARMEN and NOMAD ||^, the reactor z/g dis¬ 
appearance experiments Bugey [^ and CHOOZ and the disappearance experiment 

CDHS 1^. We will perform a fit to these data for (3-1-1) and (2-1-2) mass spectra in a unified 
formalism, which allows to compare directly the quality of the fit for these rather different 
mass schemes. 

The plan of the paper is as follows. In Sec. H we define our notation. In Secs. BE and 
|V] we consider the mixing parameters relevant in the different classes of experiments (SBL, 
solar and atmospheric, respectively), discuss constraints on these parameters and describe 
the experimental data used in the analysis. In Sec. ^ we give a thorough discussion of the 
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Figure 1: The six types of four-neutrino mass spectra. The different distances between the masses 
on the vertical axes symbolize the different scales of mass-squared differences required to explain 
solar, atmospheric and LSND data with neutrino oscillations. 


parametrization of four-neutrino mass schemes. The parameterization we introduce is based 
on physically relevant quantities and is convenient for the combined analysis. In Sec. 
we compare the (3+1) and (2+2) mass schemes when data from solar and atmospheric 
neutrino experiments are combined, whereas in Sec. |VII1| we consider the combination of 
atmospheric and SBL experiments. In Sec. ^ we present the main result of this work: a £t 
of the global set of neutrino oscillation data (solar, atmospheric and SBL) in the framework 
of four-neutrino mass spectra. We also comment on a Standard Model £t, i.e. a £t for the 
(3+0) case. In Secs. IvIIM we focus mainly on the relative comparison of (3+1) and (2+2) 
mass schemes; we will make some comments on the absolute quality of the £t in Sec. 
Finally, we present our conclusions in Sec. 


For readers interested mainly in the results of our work we suggest to skip Secs. 0- 
0 After having a look at Fig. |^, where the parameter structure of the four-neutrino £t is 
illustrated, we recommend to proceed directly to Sec. |VI|. 


II. FOUR-NEUTRINO OSCILLATION PARAMETERS 


To obtain a four-neutrino scenario from a gauge model of the weak interaction one needs to 
extend the lepton sector by a number m of 5(7(2) 0(7(1) singlet leptons |^. In such scheme 
the charged current leptonic weak interaction is specihed as by a rectangular 3 x (3 + m) 
lepton mixing matrix K = QU which comes from diagonalizing separately the 3x3 charged 
lepton mass matrix (via O) as well as the, in general (3+m) x (3+m) Majorana, neutrino mass 
matrix (via (7). Moreover the weak neutral current couplings of mass-eigenstate neutrinos is 
characterized by a non-trivial (3 + m) x (3 + m) coupling matrix P = whose effects 

will not be relevant for us as neutrinos are both produced and detected through charged 
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current interactions. If these extra singlets are all super-heavy (one example is the standard 
seesaw scheme, where m = 3), they decouple, leaving to a nearly unitary 3x3 lepton mixing 
matrix K while the projective (3 -|- m) x (3 -f m) matrix P becomes approximately the 3x3 
unit matrix (approximate GIM mechanism). 


From here on-wards we assume that, due to some symmetry or another reason [|13i |TJ, 


T5| , [T^, ^ one of the SU{2) 0 f/(l) singlets remains light enough so that it can take part 
in the oscillation phenomenology and thereby account for the LSND data. The minimum 
possibility is to have just one such light singlet, m = 1, called sterile neutrino. 

In general the physics of four-neutrino oscillations involves 3 mass-squared differences 
and the elements of the mixing matrix K. The latter have been characterized in a model- 


independent way in where an explicit parametrization was given which is, up to factor 
ordering, the standard one. In full generality K contains 6 mixing angles and 3 physical 
phases which could lead to CP violation in the oscillation phenomena [^. For convenience 
this 3x4 matrix K connecting the 4 neutrino mass helds z/j and the 3 flavor helds can be 
completed with an extra line (relating the sterile neutrino Ug to the mass eigenstates) so to 
obtain a 4 x 4 unitary matrix. In a basis where the charged lepton mass matrix is diagonal 
this leads to the matrix U diagonalizing the neutrino mass matrix: 


r'o, ^ ^ Uai^i 
i=l 


{a = e, /i, r, s). 


( 1 ) 


Because of the strong hierarchy of the mass-squared differences required by the experi¬ 
mental data the CP-violating effects are expected to be small in the experiments we consider. 
However, CP violation can be important in four-neutrino schemes for future long-baseline 
experiments, such as neutrino factories [^. Thus, neglecting the complex phases we are 
left altogether with nine parameters relevant for the description of CP conserving neutrino 
oscillations in a four-neutrino scheme: 6 mixing angles contained in U and 3 mass-squared 
differences. In the following sections we will present a choice for these parameters, which 
is convenient for the combined analysis of the different experiments and which is motivated 
by their physical interpretation. 

We label the neutrino masses as indicated in Fig. ^ and dehne for all schemes^ 


^"^LSND = m^-m^ 


and 


^”^ATM = mg - > 0 . 


( 2 ) 


All experiments we consider are insensitive to the sign of This implies that the 

(3-|-l)(j scheme is equivalent to (3-1-1)^, while (3-|-l)f, is equivalent to (3-|-l)c and {2+2)a is 
equivalent to {2+2)b-‘^ Hence, without loss of generality we can restrict ourselves to the 


^ Note that our labeling is different from the one in previous publications |^, ^ |3^. However, this 

way of labeling neutrino masses is particularly convenient as it enables a combined treatment of all the 
schemes in the same footing. 

^ These degeneracies can be lifted by considering the effects in tritium /3-decay experiments Q or 
neutrino-less double /3-decay experiments [Q. 
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discussion of the schemes (3+l)a, (3+l)f, and (2+2)^, and we always have Am^sNo > 0- The 
structure of the neutrino mass eigenstates z/ 2 , r's and z /4 is common for all these schemes. 
Only the “solar mass state” vi is inserted in different places. Let us define the index © for 
the different schemes as 


(3+l)„:0 = 2, (3+1), : 0 = 3, (2+2) : © = 4. (3) 


Then the solar mass-splitting can be written for all schemes as 


Am^oL = > 0 . 


( 4 ) 


One advantage of the labeling introduced above is that we can use the parameter Amh = 
m‘1 — mf to relate the different schemes in a continuous way. The values of Aml^ which 
correspond to the three schemes are given by 

(3+l)a : Am^i = Am^s^n + , 

(3+l)fe : Am^i = + Aml^^ - Aml^^ , (5) 

(2+2) A : Aml^ = Arngo^ • 


It will be useful to factorize the mixing matrix U into two matrices: U = 
Neglecting the complex phases in U we write the matrices 0^4 ^ product of rotation 

matrices Rij in the (i,j) subspace with the angle 9ij. We dehne 


0(1) 


^ C 14 C 13 C 12 

Cl4Cl3'Sl2 

C14 5i3 

5i4^ 

-512 

+2 

0 

0 

— 513+2 

— 5i35i2 

+3 

0 

\— 514 C 13 C 12 

— 5i4Ci3'S12 

—5i45i3 

+ 4 / 


( 6 ) 


0(2) 


R34R24R23 — 


/1 0 

0 C 24 C 23 

0 —<§34 524 C23 ~ C34523 
\0 —C34524C23 + 534 S23 


0 0 \ 

C24523 524 

— 534524523 + C34C23 534 C24 

— C34524523 — 534 C23 C34C24/ 


( 7 ) 


and order the flavor eigenstates in such a way that if all angles are zero we have the 
correspondence (i2e,i2^,i2T,i2s) = (z^i, ^’ 2 , ^^ 3 , ^^ 4 )- In the following sections we will consider 
the mixing parameters relevant in the three different classes of experiments (SBL, solar and 
atmospheric) in more detail. 


III. SBL EXPERIMENTS 
A. SBL parameters 

In SBL experiments it is a good approximation to set the solar and atmospheric mass- 
splittings to zero. Obviously, under this assumption the two schemes (3+1)^ and (3+1), 
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become equivalent. Let us define the parameters da (a = e,^,T,s) and for the two 
schemes as 


(3+1): da — lUail"^ , ^/i;e — 4 |f/e4nfo/i4p 5 

(2+2) : da = \Uai\^ + |f/„4r, = 4 \UeiU;, + 


(-) 


Then for both schemes the probability of SBL u 


(-) 


Up transitions relevant for the acceler¬ 


ator experiments LSND, KARMEN and NOMAD is given by 


Pu 


= Pu 


.Pe = sm 


2 

AE 


( 9 ) 


and the survival probabilities relevant in the SBL disappearance experiments Bugey and 
CDHS are given by 




= Pp_ 


A m2 T 

.„„ = l-4d„(l-d„)3m" 


( 10 ) 


where a = e refers to the Bugey and a = to the CDHS experiment. Here L is the distance 
between source and detector and E is the neutrino energy. It is straightforward to see that 
in the (3+1) scheme the relation 


(3 + 1): A^.e = Aded^ (11) 

holds. Hence, there are only two independent SBL mixing parameters in this case. However, 
in the (2+2) scheme the situation is qualitatively different and there is only the restriction 

(2+2): < 4min[4d^, (1-p)(l-d^)] (12) 


which follows from unitarity of U, and therefore there remain three independent mixing 
parameters for SBL experiments in the (2+2) scheme. 

Note that the probabilities Eqs. (^) and ([To|) have the same form as in the two-neutrino 
case [|I^ The amplitude A^;e can therefore be identified with the LSND mixing angle: 

A^-e = sin 2dLSND ) (13) 


and for the disappearance parameters the identification 4de(l — dp) sin^26'Bugey (and 
similar for CDHS) can be made. 


B. Constraints form Bugey and CHOOZ 

Let us consider the constraints from reactor Pg disappearance experiments Bugey and 
CHOOZ.^ To this purpose we introduce the parameter 

r]e=\Uel\^+\Upe\\ (14) 

^ Note that the Palo Verde reactor experiment Q| obtains a bound comparable to CHOOZ. As the exact 
value of this bound has very little impact on our analysis we include for simplicity only the result of 
CHOOZ. 
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which describes the fraction of the electron neutrino in the “solar sector” and is related to 
de by 


(3+1) : rje + de ^ 1, 

( 2 + 2 ) : Tje = de ■ 


The requirement that the electron neutrino must participate in oscillations with Am^QL 
in order to explain the solar neutrino anomaly leads io 7]e ~ 1. The result of the Bugey 
experiment |^6|] constrains the combination 4 de{\- — de) to be very small. Taking into account 
Eq. (p+D and //e ~ 1 one obtains |^, ^ 


(3+1) : 

( 2 + 2 ) : 


de 1 

I - de) 


< 2 X 10-2 


at 90% C.L. 


(16) 


in the relevant range of 

The disappearance probability in the CHOOZ experiment^ can be written as 


Am^ L 

^cHooz = 1-2 4(1 - de) - Achooz sin2 —, (17) 

with 

(3+1) : ^OHOOZ = 4 ? 7 e(l — de — Tje ) , flSi 

(2+2): AcHOOZ=4|f/e2r|17e3r- 

We use the result of this experiment in two ways. First, we constrain the SBL parameter 4 
similar to Bugey as described in Ref. . Second, also the parameter Achooz is constrained 
to small values. Comparing Eqs. (0) and ( |T^) and noting that for the (2+2) schemes 
(1 — 4) = |h^e 2 p + l^esP one can see that Achooz is very small in this case because of 
the bound on (1 — de ) implied by Bugey. However, for the (3+1) schemes the additional 
information of CHOOZ is important. Taking again into account that the electron neutrino 
must have signihcant mixing with z/i and z/© to obtain solar neutrino oscillations {rje ~ 1), 
we obtain the bound 


(3+1) : 1 - 4 - < 4 X 10-2 at 90% C.L. (19) 


for the values of preferred by atmospheric neutrino experiments. 

To summarize, if oscillations of Ve with Amg^^^ are required, bounds from reactor experi¬ 
ments imply in both types of mass schemes that rje has to be close to 1: for (2+2) Eqs. (0 
and (|T^ imply that (1 — rje) is bounded by Bugey, whereas for (3+1) we obtain a somewhat 
weaker bound resulting from a combination of the bounds from Bugey Eq. (p+f) and CHOOZ 
Eq. (0): 


(3+1) : 1 - r/e < 6 X 10-2 , 

( 2 + 2 ) : l-rje<2x 10-2 . 


( 20 ) 


^ Due to its value of {E/L) CHOOZ is sensitive to oscillations with rather than with AttIlskdi 

therefore it is considered as a long-baseline experiment. However, it turns out to be convenient to treat 
it together with the SBL experiments in the analysis. 
















These bounds can also be translated into bounds on the mixing angles contained in the 
factor of the leptonic mixing matrix. One of three angles in the matrix is the solar 
angle which has to be large in order to account for the results of solar neutrino 
oscillation experiments 0,0, s i i- Then the bounds given in Eq. (pOD imply in all mass 
schemes that the other two angles have to be small. 


C. Data used from SBL experiments 


In this section we describe the experimental data from SBL experiments which we are 
using for our statistical analysis. We divide the y^-function describing the SBL experiments 
into two parts: 

XsBl('^^LSND) ^LSND) de, 

Xnev(^^lsnd 5 ^LSND, de, d^) + AXlsnd(^^Lnd) ^lsno) • 

Here xIev contains the information from the experiments Bugey, CDHS, KARMEN, NO¬ 
MAD and CHOOZ, which find no evidence (NEV) for neutrino oscillations, while y: 


( 21 ) 


LSND 


includes the information of LSND, which is the only SBL experiment reporting an evidence 
for oscillations. For what concerns the parameter dependence shown in Eq. (1^) one has to 
keep in mind that in the (3-1-1) scheme ^lsno is related to dg and d^ via Eqs. o and (pAl), 
but in the context of (2-1-2) schemes all these three parameters are independent. 

The Bugey experiment searches for z/g disappearance at the distances 15 m, 40 m 
and 95 m away from a nuclear reactor. As input data for our analysis we use Fig. 17 of 
Ref. 1^^, where the ratios of the observed events to the number of expected events in case of 
no oscillations are shown in 25 bins in positron energy for the positions 15 m and 40 m, and 
10 bins for the position 95 m. The CDHS experiment [^| searches for disappearance by 
comparing the number of events in the so-called back and front detectors at the distances 
Tback = 885 m and Lfront = 130 m, respectively, from the neutrino source. The data is given 


in Tab. 1 of Ref. as ratios of these event numbers in 15 bins of “projected range in iron”. 
The KARMEN experiment [^] looks for z/g appearance in a beam. We use the number 
of positron events in 9 bins of positron energy as given in Fig. 2(b) of the second reference 
p4| . Our re-analysis of the experiments Bugey, CDHS and KARMEN is described in 
To include the results on the z/„ —z/g appearance channel obtained by 


m 


detail in Ref. 

the NOMAD experiment [Q we perform an analysis similar to the one of KARMEN. We 
use the 14 data points of the energy spectrum of charged current events given in Fig. 2 
of the hrst reference in . 


We include the result of the CHOOZ experiment ||^ by means of the y^-function 


Tchooz 


((-fcHOOz) Pexp)^ 


+ c'? 


( 22 ) 


syst 


where Pgxp = 1-01, Ustat = 2.8%, cxsyst = 2.7% [|3] and Pchooz is given in Eq. (|3). In 
the (2-1-2) case we adopt the approximation Achooz = 0 (see Sec. |HI B|) and hence y; 


2 

CHOOZ 
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depends only on the parameter de- For the (3+1) case Xchooz depends on the two independent 
parameters de and r/e- Apart from the requirement r/e ~ 1 we are not interested in the exact 
value of this parameter and we will always minimize with respect to it. In our approximation 
the CHOOZ experiment is the only one sensitive to the small value (1 — rje) (see following 
sections), and therefore the minimization with respect to rje is trivial and yields Achooz = 0. 
Again we are using only the information on de from CHOOZ in Eq. (^), which is independent 
of Therefore, the dependence on rje is not shown in Eq. (|^) . 

The total number of data points for all NEV experiments is 


A^mrv = 60 


(Bugey) 


+ 15 


(CDHS) 


+ 9 


(KARMEN) 


+ 14 


(NOMAD) 


+ 1 


(CHOOZ) 


= 99 


(23) 


To include the detailed structure of the LSND experiment [^, the LSND collabora¬ 


tion [^] has provided us with a table of the likelihood function obtained in the hnal analysis 
of their data [Hj as a function of the two-neutrino parameters ArriLSND s-nd sin^ 20 lsnd- Con¬ 
tours of this likelihood function corresponding to 90% and 99% C.L. are shown in Fig. 27 
of Ref. The reason that we can use this likelihood function, which was obtained in a 


two-neutrino analysis, also in the four-neutrino case is that the relevant four-neutrino prob¬ 
ability Eq. (^ has the same form as the two-neutrino probability. We include the LSND 
likelihood function in our analysis by transforming it into a y^-function according to 
= const — 2 ln£. Because of the event-by-event based likelihood analysis performed by 
the LSND collaboration we cannot use any information on the absolute value of the 
function. Therefore, as indicated already in Eq. (^T]) , we use in our analysis only the 
relative to its minimum: 


AXlSNd(A’TIlsnD) ^LSNd) — 2 InDj^g^p 2 InDLSND(AmLSND) ^LSNd) • (24) 

In this way we are able to include the LSND data in an optimal way, as we are using directly 
the analysis performed by the experimental group. 


IV. SOLAR NEUTRINO EXPERIMENTS 


For solar neutrino oscillations it is a good approximation to work in the limit AruLg 


CX) 


and cx), so that oscillations induced by the LSND and atmospheric mass-squared 

solar neutrino 


differences are completely averaged out. Moreover, in Refs. ^ 


oscillations in (2+2) schemes has been studied using the approximation rje = 1, which is 


justified by the Bugey bound Eq. ([Tq). The results obtained there can be applied also 
to (3+1) mass schemes, if again rje = I is adopted. Note however, that in this case only 
the somewhat weaker bound shown in Eq. (|^) applies. The solar oscillation probabilities 
obtained in these works are valid up to terms of order (1 — deY for (2+2) and [dl, (1 —rje)'^] 
for (3+1). Setting rje = 1 reduces the matrix in all cases to Riq, eliminating the other 
two mixing angles. 
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Under these approximations solar neutrino oscillations do not distinguish between (3+1) 


and ( 2 + 2 ) schemes, and depend only on the three parameters 0 sol and rj, 

solar mixing angle 6 *sol = 6^10 is given by 

tan 6 am, = 


\Uee\^ 

\U. 


The 


(25) 


ell 


and corresponds to 612 in the notation of Refs. |32 


it can be taken in the interval 


0 < 6 sol I without loss of generality. The parameter r]s is dehned by 

Vs = + |R^s0p 


(26) 


and corresponds to in the notation of Refs. 32, 33]. This parameter describes the 


fraction of the sterile neutrino participating in solar neutrino oscillations: for = 0 solar 
electron neutrinos oscillate only into active neutrinos, whereas t]s = 1 corresponds to pure 
Ue —>■ Ug oscillations.® Thus this mixing-type parameter can be interpreted as a model 
parameter interpolating between the approximate forms for the leptonic mixing matrix given 
in Refs. [jT^ and respectively. 

To include the information from solar neutrino experiments in our analysis we make use of 
the results obtained in the four-neutrino analysis performed in Ref. [Q. The experimental 
data used in this work is the solar neutrino rate of the chlorine experiment Homestake [Q, 
the weighted average rate of the gallium experiments SAGE [§], GALLEX and GNO [|], as 
well as the 1258-day Super-Kamiokande data sample in form of the recoil electron energy 
spectrum for both day and night periods, each of them given in 19 data bins, and the recent 
result from the charged current event rate at SNO ^j. The total number of data points 
contained in xIol 


Afam, — 3 


(Cl,Ga,SNO) 


+ 38 


(SK) 


= 41 


(27) 


Details of the solar neutrino analysis can be found in Refs. ||31| , |33| , pj and references therein. 

To include the results of Ref. in our analysis we use Xsoh ^ function of rjs (minimized 
with respect to the other two parameters and 6 'sol) shown in Fig. 3 of Ref. |^, which 

we reproduce in Fig. |^. The is shown relative to the global minimum, which lies in the 
large mixing angle (LMA) region and has the value (XsoL)min = 35.3 for Nsql — 3 = 38 
degrees of freedom (d.o.f.). The three lines in the hgure are obtained by requiring that 
the solution of the solar neutrino problem lies in the three regions LMA, low/quasi-vacuum 
(LOW) and small mixing angle (SMA), respectively. Note that XsohiVs) is the same for all 
mass schemes. We clearly see from this hgure that solar neutrino data prefers rjg = 0, i.e. 


pure active oscillations. At 99% G.L. there is the upper bound from the solar data [33]: 


solar data: ps < 0.52 . 


( 28 ) 


® The parameter rjs is similar to the parameters + and c^, which have been introduced in Refs, and 
[p0[, respectively, to describe the effect of (2+2) mass schemes in Big-Bang nucleosynthesis. 
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Figure 2: 
presented 


Axsol cis a function of rjs for the different solutions to the solar neutrino problem, as 


in Fig. 3 of Ref. 133(1 . 


V. ATMOSPHERIC NEUTRINO EXPERIMENTS 


For the oscillations of atmospheric neutrinos it is a good approximation to set Arrig^^ to 


zero and to also assume the limit Am?„ 


oo. 


In Refs. 


52| hts of atmospheric 


neutrino data in a (2+2) framework have been performed by making use of the Bugey 
constraint Eq. ( p+Q and setting rje = 1. The approximations ArngQ^ = 0 and rje = 1 imply 
that the electron neutrino decouples completely from atmospheric neutrino oscillations. In 
(3+1) spectra the contribution of electron neutrinos to atmospheric oscillations is limited 
by the somewhat weaker bound shown in Eq. (p0|) ; however, in Ref. it was found that 
a z/g contamination small enough not to spoil the result of the CHOOZ experiment has only 
a very small effect on the quality of the £t of atmospheric neutrino data. Therefore, it is 
justihed to adopt the approximation r]f. = 1 also for (3+1) schemes . 


Under these assumptions, atmospheric neutrino oscillations do not distinguish between 
(3+1) and (2+2) schemes, and reduce to an effective three-neutrino problem involving only 
the flavors z/^, z/^, the mass eigenstates z/ 2 , r's, and the mixing matrix dehned i 


m 


Eq. (]^). The function depends on the four parameters 623 , 6*34 and 6*24 |^, |55| , 

and to cover the full physical parameter space one can choose the ranges 0 < (6*24, 6*34) < 7 r /2 
and — 71/2 < 623 < vr/2.® Therefore, in addition to the two parameters Arril^^ and ^atm = 
023 corresponding to the two-neutrino parameters, we need two more angles to describe 


® Note that ( 034 , 6 * 24 ,^ 23 ) in our notation correspond to (z? 24 , 023, ^ 34 ), respectively, in the notation of 
Refs. 
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atmospheric neutrino oscillations in a four-neutrino framework . 

To understand the physical meaning of the angles 6^24 and 634 let us consider their relation 
to the parameters and dg, which we have dehned in Eq. (H). Under the approximation 
r^e = 1 , we obtain in all the schemes 


^ _ |(o( 2)|2 _ 2 

^ _ 1 ^( 2 ) |2 _ 2 2 
Ug — \'^g4 \ — <- 24'-34 • 


(29) 


The quantity (1 — d^) [(1 — d^)] corresponds to the fraction of the muon [sterile] neutrino 
participating in “atmospheric” neutrino oscillations. For d^ = S 24 = 0 the muon neutrino 
lies completely in the atmospheric sector, while for the (strongly disfavored) case d^ = 1 
there are no oscillations of with the scale Hence, atmospheric data will constrain 

dn to be small. Depending on the value of the bound on dn is strengthened by the 


(-) 


SBL disappearance experiment CDHS 


Similarly, dg = = 1 corresponds to 


pure active atmospheric oscillations, whereas for d^ = 0 the sterile neutrino fully participates 
in oscillations with The cases correspond to the approximations used in the early 


m- 


papers Big: the mixing-type parameter dg can be interpreted as a model parameter 
terpolating between the approximate forms for the leptonic mixing matrix given in Refs. 

{dg = 0 ) and {dg = 1 ), respectively. 

For the atmospheric data analysis we use the following data from the Super-Kamiokande 
experiment |0|: e-like and /i-like data samples of sub- and multi-GeV, each given as a hve- 
bin zenith-angle distribution, up-going muon data including the stopping (5 bins in zenith 
angle) and through-going (10 angular bins) muon fluxes. Further, we use the recent update 
of the MACRO [Q up-going muon sample (10 angular bins). We obtain a total number of 
data points contained in x\tm of 


Aatm — 35 (SK) + 10 (MACRO) — 45 . 


(30) 


For further details of the atmospheric neutrino analysis see Refs. [^, |^, ^ and references 
therein. 

In Fig. 1^ we show the results of our atmospheric neutrino analysis regarding the angles 
024 and 6 * 34 . This hgure corresponds to Fig. 6 of Ref. |^, but now using the updated results 
of MACRO. In the upper panel we show the 90% and 99% C.L. allowed regions (2 d.o.f.) 
for the parameters dg = C 24 C 34 and d^ = S 24 . To obtain these regions we minimize x\tm 
with respect to the other two parameters ^atm and As expected, atmospheric data 

constrains d^ to small values, implying a large fraction of 12 ^ participating in atmospheric 
oscillations. For what concerns the parameter dg, values close to 1 are preferred, which 
means that 12 ^ oscillates mainly to active neutrinos. This can be seen clearly from the lower 
panel of Fig. where we display Ax2^,^(4) = xItm^^s) - {xItm)^:^ ■ Here (xiTM)min = 27.9 
for Aatm — 4 = 41 d.o.f. and Xatm is minimized with respect to all other parameters. We 
show the line for the updated MACRO data P| and compare with the line obtained from 
the old MACRO data [0] , which corresponds to the data used in Ref. . For large values 
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Figure 3; (a) 90% and 99% C.L. allowed regions for the parameters dg = C 24 C 34 (ordinate) and 
d^ = S 24 (abscissa) from atmospheric neutrino data. The best fit point is marked with a star, (b) 
Axatm S’® a function of dg using old and new MACRO data. Also shown are the Ax^-values 
corresponding to 90% and 99% C.L. for 1 d.o.f.. 


data set 

parameters 

solar 

> ^SOL , f]s 

atmospheric 

Am^TM ) ^ATM) ^24 ) ^34 

SBL appearance 

) ^LSND 

SBL disappearance 

^^LSND 5 


Table I: Four-neutrino parameters for the different data sets. 

of dg the lines are very similar, however for small values the £t gets worse. This means 
that atmospheric data get stronger in rejecting a sterile component in atmospheric neutrino 
oscillations. 

VI. FOUR-NEUTRINO PARAMETERS IN THE COMBINED ANALYSIS 

In the previous sections we have discussed the parameterization of the four-neutrino 
problem for the different data sets separately. We summarize our choice of parameters in 
Tab. |. Note that for (3-1-1) schemes rje is an additional independent parameter, but we do 
not list it in Tab. | because in our approximation CHOOZ is the only experiment sensitive 
to it and we always minimize with respect to it. In (2-1-2) schemes we have r/e = d^ according 
to Eq. (1^. We have chosen the parameters listed in Tab. | in such a way that they have 
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Figure 4: Parameter dependence of the three data sets solar, atmospheric and SBL. Exact defini¬ 
tions of the parameters are given in Secs. g0and0. 


a well-defined physical meaning in the context of a given data set. Note that this physical 
interpretation is independent of the mass scheme: for example, regardless of whether we 
assume (3-1-1) or (2-1-2) schemes, ps is the fraction of sterile neutrinos in solar oscillations, 
6*24 describes the fraction of in atmospheric oscillations (see Eq. (p9|)), and sin^26*LSND is 
the SBL —>■ Ve amplitude, and so on. The fact that it is possible to describe the results of 
any of the given set of experiments in terms of physical quantities independent of the mass 
scheme implies that none of the considered data sets (solar, atmospheric, SBL appearance 
or SBL disappearance) can be used on its own to distinguish between different mass spectra. 
This follows from the approximation r/e ~ 1, which is motivated by the bounds from reactor 
neutrino experiments, and from the strong hierarchy among the mass-squared differences 
indicated by the data. This hierarchy implies that for any set of experiments only one mass 
scale is relevant. In the following we show in detail that the differences between the mass 
schemes manifest themselves only if two or more data sets are combined, i.e. if the relation 
among parameters belonging to different data sets is considered. 

For the combined analysis we will describe neutrino oscillations by means of the following 
parameters: beside the three mass-squared differences ArngQ^; ^^atm 'w® use 

the six parameters 

^SOL 5 ^ATM ) ^LSND 5 Vs 5 5 Ve ■ ( 31 ) 


It is easy to check that indeed for all mass schemes these six parameters - dehned as in the 
previous sections - can be used to describe, in a physically more convenient way, the most 
general CP conserving leptonic mixing matrix |^. Each of the y^-functions describing the 
three data sets (SBL, solar, atmospheric) depends only on a sub-set of these parameters: 


TsOl('^^SOL) ^SOL) hs) ) XaTm(^^ATM) ^ATM) d^, Tjg) , ]yggj^(A?7lj^gj„p, ^lSND) d^, T]g) . (32) 


We illustrate the parameter dependence of the data sets in Fig. |. The three angles 6 *sol 5 
6*ATM, 6*lsnd are related directly to the amplitude of the oscillations in the corresponding 
experiments solar, atmospheric and LSND, while the two quantities rjs and d^ account for 
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the coupling between different data sets. As indicated, the parameter r]s is in common to 
solar and atmospheric neutrino oscillations; if we express it in terms of the atmospheric 
angles we obtain the different relations, depending on the mass schemes: 

(3+l)a • Vs — ~ ('S 24 C 34 CATM 'S34Satm) j 

(3 + l)b : Vs = ~ ('S24C34’Satm + ’S34Catm)^ , (33) 

(2+2): r,. = |02P = 44. 


On the other hand the parameter is in common to SBL and atmospheric neutrino oscilla¬ 
tions; here the coupling is the same in all mass schemes (see Eq. (0))- Another important 
difference between (3-1-1) and (2-1-2) arises due to the combination of SBL appearance and 
disappearance experiments (see Eqs. (0) and ([T^) ). There is no direct coupling between so¬ 
lar and SBL oscillations, they do not depend on a common parameter. This simple coupling 
of the data sets is a nice feature of our parameterization (within the adopted approxima¬ 
tions) which renders the combined analysis possible, despite the large number of parameters 
involved. Note that only the SBL experiments involve the additional parameter rje explicitly, 
since - as stated in Secs. 0 and 0 - for what concerns the analysis of the other two data 
sets it is safe to assume r^e = 1. In the (2-1-2) scheme we have the relation rje = de and this 
parameter is important for the SBL disappearance amplitude, whereas in (3-1-1) only the 
long baseline reactor experiment CHOOZ is sensitive to Ve and we always minimize with 
respect to it. 

For the unified analysis of all the mass schemes we will consider and 

Aml-^ as the three independent mass-squared differences. The as a function of 
will have three local minima corresponding to the schemes (3-|-l)a, (3-|-l)f, and (2-\-2)a at 
the values given in Eq. (^. Beside this parameter indicating the scheme we will display the 
results of our numerical analysis using the following parameters. For the analysis of solar 
and atmospheric data in Sec. [VII| we consider the as a function of Vs', the results of the 
analysis of atmospheric and SBL data (Sec. [VII1| ) are given in the (Am^sNO) 26 'lsnd) plane, 
while for the fully global analysis (Sec. m we use all three parameters (AruLSN^, 6 *lsnd, hs)- 
Note that in a x^ analysis the size of the allowed regions depend crucially on the number of 
parameters considered. Our aim was to identify which parameters describe the most relevant 
features of the physics in each case. 

Before closing this section let us note some subtleties related to our parameterization. 
As described above, some of the parameters shown in Eq. (^Ij), which we are using for our 
global fit, obey different relations depending on the mass scheme considered. The question 
arises of how to treat these different relations among parameters in a common framework 
for the two mass schemes. Indeed, we are using the parameter Am‘li to formally describe a 
continuous transition between the vastly different mass spectra. 

Let us consider a completely arbitrary parameterization of the general four-neutrino prob¬ 
lem p9|. We have 3 mass-squared differences and 6 angles, e.g. the angles 9ij introduced in 
Eqs. (^ and (0). Now one can think of a fit to the data in this general parameterization. 
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Figure 5: Ax^tm a function of the fraction of the sterile neutrino in solar oscillations rjs for all 
four-neutrino mass schemes, (3-|-l)a, (3-|-l)f, and (2-1-2). 

In practice it is almost impossible to perform this general nine-parameter fit with current 
computer technology. However, the results of such an analysis would be six well separated 
regions in the nine-dimensional parameter space, corresponding to the six mass schemes 
shown in Fig. |^. In these islands our approximations hold and it makes sense to adopt a 
parameterization motivated by phenomenology. The allowed regions for the parameters 9ij 
can be mapped to allowed regions for the parameters shown in Eq. (^), which have a simple 
physical interpretation. Inside any given island it is clear which relations among the new 
parameters have to be applied. Obviously this reasoning is only valid under the assumption 
that various regions are well separated. This assumption can be justihed by noting that 
to move continuously from a (2-1-2) to a (3-1-1) scheme one has to break up the hierarchy 
among the mass-squared differences, and we can expect that at least one data set will give 
a large which separates the corresponding allowed regions. 


VII. ANALYSIS OF SOLAR AND ATMOSPHERIC NEUTRINO DATA 


In this section we combine solar and atmospheric neutrino data. In Ref. these data 
have been considered in the (2-1-2) scheme. Here we discuss some slight changes in this case 
(due to the updated MACRO data), we extend the analysis also to the case of the (3-1-1) 
mass scheme in a way that allows a direct comparison of the £t for these two schemes. 

Before combining the two data sets let us consider the impact of atmospheric data alone 
on the parameter t]s, describing the fraction of sterile neutrinos in solar oscillations. The 
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relation of rjg to the atmospheric parameters is given in Eq. for the three mass schemes 
(3+l)a, (3+1)6 and (2+2). In Fig. | we show AxiTM(h^) = xItm^Vs) - (xiTM)min for the 
three cases, minimizing with respect to the other parameters 6 *atm and d^. The line 

corresponding to the (2+2) scheme is identical to the one shown in the lower panel of Fig. 
because of the ( 2 + 2 ) relation rjg = dg = C 24 C 34 (see Eq. (0)). Atmospheric data prefers 
large values of t]s, which corresponds to active atmospheric neutrino oscillations. From 
the hgure we can read off the 99% C.L. bound 


atmospheric data: rjs > 0.54 for (2+2) schemes 


(34) 


which is in disagreement with the bound (^) from solar data. On the other hand, concerning 
(3+1) schemes in Ref. the very interesting fact was noted that atmospheric data give 
a constraint on the fraction of the sterile neutrino participating in solar oscillations. From 
Eq. (1^) it follows that = 1 — is the fraction of sterile neutrinos in 

atmospheric oscillations which should be small according to the data |^. Comparing this 
with Eq. (^3l) we expect that for (3+1) schemes atmospheric data prefers small values of rjg. 
Indeed, from Fig. ||we hnd the 99% C.L. bounds 


atmospheric data: 


r]s < 0.35 
rjs < 0.42 


for (3+l)a schemes, 
for (3+l)fe schemes 


(36) 


which are even stronger than the one from solar data Eq. ([^). 

From Fig. ^ one can see that, although there are quantitative differences between the 
two schemes (3+1)^ and (3+1)^, the qualitative behavior is similar. Looking at Eq. (|33[) 
it is easy to see that the relation between rjs and the atmospheric angles in the (3+l)a 
scheme reduces to the one in the (3+1)^ scheme under the transformation ^atm —^ ^atm — 
7 r/ 2 . Such a transformation, when applied to the atmospheric oscillation probabilities, is 
equivalent to changing the sign of Am^.j,^, hence we have ^atm — 7 r/ 2 , 6 ^ 24 , 6 * 34 ) = 

Xatm(~'^^atm) ^atm, ^ 24 , 6 * 34 ). Therefore, the origin of the difference between the two schemes 
(3+l)a and (3+1)^ can be explained in two different ways. If we require > 0, then 

we end up with two different relations between rjg and the atmospheric angles ^atm, ^24 and 
034 , as we have done so far. Alternatively, if the dehnition of the atmospheric angles is 
adjusted so that their relation with rjg is the same in (3+1)^ and (3+1)^ schemes, then it is 
no longer sufficient to restrict to the case > 0 , and the case < 0 should be 

investigated as well. In the latter approach, it is clear that the difference arise because of 
the presence of matter effects in atmospheric neutrino oscillations, which are sensitive to the 
sign of Since in this work we are mainly interested in the comparison of ( 2 + 2 ) with 

(3+1) in general, from now on we will always minimize with respect to (3+1)^ and (3+1)6 
by choosing from the two corresponding values of Am^^ the one with the lower 

From Eqs. (PB[) , (^41) and (|45| ) one expects that combined solar and atmospheric neutrino 
data will prefer (3+1) mass schemes over (2+2). In order to quantify this statement let us 
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Figure 6: Combined x^-function for solar and atmospheric neutrino data for (3+1) and (2+2) for 
different solar neutrino solutions with respect to the global minimum (see text for details). 

consider the following x^-function: 

XsOL+ATm(^S) ^^4l) = XsOL d" XaTM (36) 

where we minimize with respect to the parameters 9sol, 9atm and d^. As 

explained before (see Sec. II) the parameter Am\^ relates the different schemes. In Fig. ^ 
we show the Ax^ projected onto the T^s-axis for the two regions of Am'h corresponding 
to the schemes (3+1) and (2+2) according to Eq. (|]). In both cases, (3+1) and (2+2), 
we refer to the same minimum, which occurs for the (3+1) scheme with rjs = 0. For the 
dashed (dashed-dotted) line we restrict the solar solution to be LMA (LOW), while for 
the solid line (labeled “unc” in the figure) we choose the solution which gives the weakest 
restriction (unconstrained). This corresponds to our current knowledge of the solution to 
the solar neutrino problem. The reason why the line for LOW sometimes is below the one for 
the unconstrained case is that it is referred to the minimum for the LOW solution, which is 
higher than the minimum for LMA (which coincides with the minimum in the unconstrained 
case) as can be seen from Fig. For the unconstrained case the minimum has the value 
(xLL+ATM)min = 63.2 for Nsol + Aatm - 7 = 79 d.o.L. 

As expected, solar and atmospheric neutrino data prefer (3+1) because in this case both 
can be explained by active neutrino oscillations. The dotted line in Fig. ^ shows the value 
of Ax^ = 9.21 corresponding to 99% C.L. for 2 d.o.L, which are the two parameters rjs and 
Am|^. Therefore, for the unconstrained and LMA cases (2+2) is disfavored at more than 
99% C.L. with respect to (3+1). If we compare the local minimum with respect to rjs in 
(2+2) with the global minimum in (3+1) we hnd for the unconstrained case 
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Figure 7: Combination of atmospheric and SBL data. We show projections of the three-dimensional 
90% and 99% C.L. regions corresponding to (3-1-1) and (2-1-2) in the (AttIls^d; sin^ 20 lsnd) plane. 
The best fit point lies in {‘2+2) and is marked with a star, the local best fit point in (3-1-1) is marked 
with a circle. The doted line is the 99% C.L. region from LSND data alone [11|. 


Conversely, the LOW (2-1-2) solution is still allowed at the 99% C.L.. The reason for this 
is, that the LOW solution is not as strong to reject a sterile component in solar oscillations 
as the LMA solution (see Fig. so that the disagreement with atmospheric data in the 
context of (2-1-2) schemes is somewhat weaker. Let us note that for the unconstrained case 
large values of rjg ~ 0.76 (corresponding to a large component of sterile neutrino in solar 
oscillations) are slightly preferred over small ones. This means that the inclusion of the 
updated MACRO results makes atmospheric neutrino data slightly more powerful to reject 
the sterile neutrino than the unconstrained solar data. 


VIII. ANALYSIS OF ATMOSPHERIC AND SBL NEUTRINO DATA 

In this section we combine the data sets from atmospheric and SBL neutrino experiments. 
To this purpose we consider the x^-function 

TaTM-|-SBl(^^415 ^^LSND! ^LSNd) XaTM a Xsbl ■ (38) 

From Eq. (^) we can see that the terms on the right hand side in general depend on the 
parameters (Am^gNp, 6 *atm, ^'lsnd, "^^e)- To obtain the parameter dependence as 

shown in Eq. (^8]) we proceed as follows. First, we minimize x\tm with respect to Am^.j,^, 6 *atm 
and rjs. In a second step we minimize with respect to dg ^md d^ by taking into account the 
relation ([TI|) or (0). depending on the mass scheme considered (he. on the value of 
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The allowed regions in the parameter space {Amh, Am^sNO) ^lsnd) are given by = 6.3 
(11.3) for 90% (99%) C.L. (3 d.o.f.). In the right panel of Fig. ^ we show a projection of the 
three-dimensional regions corresponding to the (2-1-2) case, which include the best fit point 
(AmLSND = 0.91 eV^, sin^ 26 *lsnd = 3.16 x 10“^). One can see that the allowed regions cover a 
large part of the two-neutrino allowed region by LSND alone, which is shown as the dotted 
line. The allowed region disappears for values of sin^ 26 *lsnd ^ 0.06 because of the constraint 
from the Bugey experiment. At large values of Am^gj^n the bounds from KARMEN and 
NOMAD are important. In the left panel we show the projection of the three-dimensional 
volume corresponding to the (3-1-1) case with respect to the global minimum, which lies in 
the (2-1-2) plane. Only four small islands appear at 99% C.L.. If we compare the local best 
fit point for (3-1-1) at (Am^gp^p = 1.74 eV^, sin^ 201^3^,0 = 1.41 x 10“^) with the global best £t 
point we find 

= (xa™+sbl) - (xlimL) = 6.9. (39) 

The conclusion from Fig. ^ and Eq. (^) is that SBL data combined with atmospheric data 
clearly prefer (2-1-2) over the (3-1-1) spectra. Fig. ^is a beautiful confirmation of the results 
of our previous work |^, where we have analyzed a similar data set in the (3-1-1) framework, 
but using a very different statistical method. The reason for the (2-1-2) preference by the SBL 
data is well known |^, ^ ^ ^ and can be understood from Eqs. (P!TD and (^. 

For the (3-1-1) case both and must be small because of Bugey and CDHS, respectively, 
which leads to a double suppression of the LSND amplitude sin^ 201^3^0 according to Eq. (©• 
In contrast, for the (2-|-2) case (1 — de) and have to be small and Eq. (|l^) implies only a 
linear suppression of sin^ 2 dL 3 ND. 

IX. GLOBAL ANALYSIS 

The results of the previous section, i.e. that atmospheric-l-SBL data prefer (2-1-2) over 
(3-1-1), are in direct conflict with the results of Sec. ivg. where we have found that so¬ 
lar-(-atmospheric data prefer (3-1-1) over (2-1-2). This shows that there is some tension in 
the existing data in a four-neutrino framework, and to clarify the situation it is necessary 
to perform a combined analysis of all the data. To this end we consider the y^-function 

Xglobal(^^4l! ^^LSNDl ^LSND) Vs) XsOL “F XaTM A XsBL ' (^6) 

We minimize the right-hand side of this equation with respect to all the parameters, except 
the ones shown on the left-hand side. Allowed regions are given by 

^Xglobal = Xglobal ~ (Xglobal)min = 7.8 (13.3) (41) 

for 90% (99%) C.L. (4 d.o.f.). From this equation we obtain two four-dimensional volumes 
in the space of (Am^^, Am^gj^p, sin^ 2dL3ND, r^^) corresponding to (3-1-1) and (2-1-2), which we 
display in Fig. | in the following way. In the lower panels we show projections of the 
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Figure 8: Global combination of current neutrino oscillation data: solar, atmospheric and SBL. 
We show the AXgio^ai ^ function of r]s (upper panels) and projections of the four-dimensional 
90% and 99% C.L. regions on the (Am^g^D) sin^ 20 lsnd) plane (lower panels) (see text for details). 
The global best fit point lies in (3-I-1) and is marked with a star, the local best fit point in (2+2) 
is marked with a circle and the local minimum in (2-|-2) is marked with a triangle. The doted line 


in the lower panels is the 99% C.L. region from LSND data alone |]T]. 


four-dimensional volumes onto the (AtUlsno) sin^ 26 *lsnd) plane. In the npper panels we show 
AXgiobai minimized with respect to all parameters except rjs. The projections of the fonr- 
dimensional 90% and 99% C.L. volnmes onto the r^^-axis are given by the intersections of 
the solid lines in the npper panels with the corresponding horizontal dotted lines. 

Let ns discuss the results of the global analysis. We hnd that the global minimnm lies in 
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^"^LSND [eV^] 

sin^ 20Lg].^n 

Vs 

djjL 

(ig 

global best fit (3+1) 

1.74 

1.41 X 10-3 

0.0 

1.98 X 10-2 

1.79 X 10-2 

best fit (2+2) 

0.87 

3.55 X 10-3 

0.93 

6.56 X 10-3 

0.99275 

local minimum (2+2) 

0.87 

3.55 X 10-3 

0.21 

1.32 X 10-2 

0.99275 


Table II: Parameter values at the best fit points in (3+1) and (2+2) and at the local minimum in 
( 2 + 2 ). 


the (3+1) scheme. This minimum is marked as a star in Fig. ^ In (2+2) there are two local 
minima: the (2+2) best £t point is marked with a circle and corresponds to large values 
of r]s, whereas the second local minimum (marked with a triangle) occurs for small rjgJ 
The values of the parameters at these minima are given in Tab. However, the difference 
between the best fit points in (3+1) and (2+2) is not very big: 



(42) 


We conclude that the schemes (3+1) and (2+2) give a comparable global £t to the data. 
This can also be seen from the fact that there are large allowed regions for both mass spectra. 
The conflicting values given in Eq. (|37|) (for solar and atmospheric data) and in Eq. (|3^ 
(for atmospheric and SBL data) cancel each other to some extent. Solar plus atmospheric 
data seem to be stronger than SBL data, therefore (3+1) is slightly preferred over (2+2) in 
the global £t to current neutrino oscillation data. 

The shape of the allowed regions in the (Am^sNO) sin^ 26 'lsnd) plane for (2+2) (lower right 
panel of Fig. |^) is simliar to the one expected from a two-neutrino analysis of SBL data alone. 
In the region 0.18 eV^ < ~ § eV^ they follow closely the two-neutrino LSND region. 


The slight shift to smaller values of is because of the constraint from KARMEN. 

Large values of sin^ 26 'lsnd ^ 0.06 are excluded by Bugey and for AruLg^D ^ 10 eV^ constraints 
from KARMEN and NOMAD are important.^ In contrast, in the (3+1) case (lower left panel 
of Fig. P) the allowed regions consist of several islands and are very different from the two- 


neutrino ones. The most prominent islands are at the values Am 

sin^ 29 


LSND 


LSND 


0.9, 1.7, 6 eV^ and 

10“^. These are the values of Am^g^D where the bounds of all NEV experiments 


have some marginal overlap with the LSND allowed region [^T], |^, |^, However, 

at 99% C.L. appears an allowed region at Am^gNo ~ 0-5 eV^, and a very marginal island at 
rsj 2.5 eV^. There is also an allowed region for large values of Am^g^D ^10 eV^. 
However, in this region there are further constraints from experiments not included in our 


^ Note that the two stars in the lower and upper panels actually correspond to the same single point in 
the four-dimensional space. The same holds for the two circles. We do not show the triangle in the lower 
right panel, because it would coincide with the circle (see Tab. 

® All the relevant SBL bounds are shown e.g. in Fig. 27 of Ref. |ll[|. A combined analysis of LSND and 
KARMEN in a two-neutrino framework has been performed in Ref. ||^. 
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analysis, which are BNL E776 


h-) 


(-) 

Vp. 


appearance and disappearance). Therefore, we do not display valnes of > 20 eV^ 


appearance) and CCFR 


d-) 


(-) 

Vp 


Discnssing the npper panels of Fig. | we note that, as already fonnd in Sec. [VII|, large 
valnes of are preferred for the (2+2) case. Comparing the shape of the in the global 
analysis with the one shown in Fig. ^ we observe that the inclusion of SBL data strengthen 
this trend to some extend. This implies a large component of the sterile neutrino in so¬ 
lar neutrino oscillations and corresponds to the LOW/quasi-vacuum solution of the solar 
neutrino problem (see Fig. But also the local minimum for smaller 77 ^ values, which 
corresponds to the LMA solar solution and implies a large sterile component in atmospheric 
oscillations, is well inside the 90% C.L. region. The difference in between the two local 
minima in (2+2) is 2.7. Moreover, the minima in are not very deep so that all values of 
ris between 0 and 1 are within the 99% C.L. region. Only values around 0.5 are excluded at 
90 % C.L.. The results shown in Fig. | were obtained by using unconstrained solar data. 
We have also performed the analysis by restricting solar data to the LMA and LOW region. 
The results are very similar to the unconstrained case. For the LMA solution we obtain 
approximately the solution corresponding to the local minimum in (2+2), which means that 
(2+2) is sightly more disfavored against (3+1), whereas for the LOW case the difference 
would become even smaller than shown in Fq. (^). 

Recently, solar neutrino data have been analyzed using a new prediction of the ®B flux 


From Tab. 3 of Ref. [Q one can see that LMA becomes relatively better than LOW. This 
would lead to an up-wards shift of the LOW line in Fig. of approximately 3 units. Con¬ 
sequently solar data becomes stronger in rejecting the sterile neutrino. Regarding the four- 
neutrino analysis, this would disfavor the (2+2) scheme slightly more against the (3+1) 
case. 

In our framework it is also possible to test the ht of the (3+0) scenario, where the solar and 
atmospheric neutrino problems are explained by oscillations among three active neutrinos 
and the explanation of LSND is left out. This would correspond to the Standard Model 
situation. We obtain this case by considering the (3+1) scheme (this hxes the parameter 
Am 4 ^) and setting the parameters df, = = rjs = 0. Then the sterile neutrino is completely 

decoupled and we are left with three active neutrinos and the mass splittings ArUg^j^ and 
Am^TM- We hnd a difference in to the best ht point of 


AV (Xg,‘77) . =19.8. 


(43) 


For 4 d.o.f. {de, Amh) this corresponds to an exclusion at more than 99.9% C.L.® We 
conclude that the data of LSND (using the result of the analysis performed by the LSND 
collaboration) plays a very signihcant role and that the global ht in a four-neutrino scenario 
is much better than in the three-active neutrino case. 


Regarding the exact value of this C.L. see also the discussion of the (3-1-0) case in the next section. 
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X. GOODNESS OF FIT 


In the previous sections we have restricted ourselves to the relative comparison of the 
ht in the different mass schemes. Here we discuss the absolute goodness of £t (GOF). A 
common way of evaluating the GOF is to consider the absolute value of the y^-function at 
the best £t point. We are aware of the fact that GOF-values obtained in this way are not 
very restrictive in a global analysis with many data points like in our case. Therefore, we 
will also consider in this section the quality of the £t in the four-neutrino schemes (3-1-1), 
(2-1-2) and for the three active neutrino case (3-1-0) for each of the different data sub-sets 
separately. 

As explained in Sec. |1II Q we are not able to use any information on the absolute value 
of xLnd from the LSND data. However, let us note that the ht for LSND is expected to 
be rather good. In Ref. the fh® fit of fh® L/E distribution to the decay-at-rest 

events of the LSND data is given for two typical values of Am^gNo as = 4.9 and 5.8 for 
8 d.o.f., which corresponds to a very good GOF of 77% and 67%, respectively. From Fig. 
one can see that for (3-1-1) as well as for (2-1-2) the best ht point lies well inside the 99% 
G.L. region of LSND. Therefore, we expect that the contribution of LSND will not worsen 
the global ht significantly (see also Tab. ra¬ 
in the following we evaluate the y^-functions for all experiments except LSND 


2 = 2 , 2,2 
Xglobal— LSND — XsOL ' XaTM ' XnEV 

at the global best ht point in (3-1-1) and the best ht point in (2-1-2): 

(3+1) : x|obai-LSND = 150.0/176 d.o.f., 
(2+2) : y|„bai-LSNo = 156.1/176 d.o.f.. 


(44) 


(45) 


The number of d.o.f. is given by (see Eqs. (|^), (|^), (|30|) ) + A^tm + A^ev = 185 

minus 9 htted parameters. Usually a ht is considered to be good if the value of the is 
approximately equal to the number of d.o.f.. 

The GOF implied by the y^-values and the corresponding number of d.o.f. given in 
Eq. (43) would be excellent for both schemes. However, one has to be careful in the inter¬ 
pretation of these numbers. This y^-test for the GOF is not a very restrictive criterion in a 
global ht of diherent experiments with a large number of data points and many parameters. 
One reason is that in such a case a given parameter is constrained often only by a small 
sub-set of the data. The rest of the data (which can contain many data points) is htted 
perfectly by this parameter (because it is insensitive to it). A discussion of this problem can 


be found in Ref. or in the context of solar neutrino analysis in Refs. 


In order to obtain more insight in the quality of the global ht we will consider the following 
quantities: 


Ay^ = y^(a) - (y^). 


(46) 
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data set 

d.o.f. 

parameters 

(3+1) (2+2)best 

(2+2)iocal 

(3+0) 

solar 

3 

^'sOL, Vs 

0.0 

10.7 

1.6 

0.0 

atmospheric 

4 

A?77.atm? ^ATM; 'Hsi 

0.0 

0.2 

11.5 

0.3 

LSND 

2 

0LSND 

3.0 

0.7 

0.7 

29.0 

NEV 

2/3 

^LSND) dg, 

8.8 

3.7 

4.1 

2.3 


Table III: /S.y^ for the different data sets of the best fit points in (3+1) and (2+2), the local 
minimum in (2+2) and for the (3+0) case (see text for details). Also shown is the number of d.o.f. 
and the corresponding parameters for each data set. 


Here is the x^-function of the data set a = SOL, atm, lsnd, nev and (xCT)min is the 
minimum of x^- This quantity can be used to test if a given point a in the parameter space 
is in agreement with the data set a. For a we will use the best £t points from the global 
analysis for (3+1) and (2+2), the local minimum in (2+2) and the point corresponding to 
the (3+0) case. This approach is similar to the method proposed in Ref. [5^. 

Let us discuss the results of this analysis, which are shown in Tab. |ITIi One can see that 
solar and atmospheric data are in perfect agreement with the global best fit point in (3+1). 
The reason is that in this case both effects are explained by active neutrino oscillations, 
which is preferred by the data. Also a Ax^ = 3 for LSND is in good agreement; the best 
£t point lies within the 90% C.L. region for the two parameters sin^26'LSND- 

However, the (3+1) best £t point gives a rather bad £t to the SBL experiments hnding no 
evidence of oscillations: a value of Ax^ = 8.8 for 2 d.o.f. {dg and d^) is ruled out at 98.7% 
C.L.. 


Regarding the (2+2) scheme we observe some problems in the fit of solar or atmospheric 
data: At the best £t (2+2) solution we obtain for solar data Ax^ = 10.7 for 3 d.o.f., 
which is ruled out at 98.7% C.L., while the fit of the other data sets ATM, LSND and 
NEV is very good. The reason for the problems in the solar data is that the best £t for 
(2+2) prefers a large value of r]s (corresponding to the LOW solution). This implies a large 
component of the sterile neutrino in solar oscillations, which gives a bad £t. In the local 
minimum for (2+2) ~ which is marked with a triangle in the upper right panel of Fig. ^ and 
corresponds to the LMA solution - the £t of solar data is very good, whereas atmospheric 
data gives a Ax^ = 11-5 for 4 d.o.f., which is ruled out at 97.8% C.L.. In this case the 
bad fit is a consequence of the large sterile component in atmospheric oscillations implied 
by the small value of rj^. The interesting shape of Xgiobai(hs) (2+2), which disfavors equal 
sterile admixture in solar and atmospheric oscillations, implies that either the solar or the 
atmospheric £t is bad in the (2+2) case, but never both. 

From the last column in Tab. H one can see that all experiments except LSND are in 
perfect agreement with the (3+0) scenario. However, as expected the fit of LSND is very 
bad in this case and yields a Ay^ = 29 for 2 d.o.f.. In the Gaussian approximation implied 
by Eq. (^^ this would be ruled out at an extremely high C.L.. Let us note that far from 
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the allowed region of LSND this approximation may not be completely justihed. However, 
it is evident that LSND data is in strong disagreement with no oscillations. According to 
Table X of Ref. Il^ the probability that the observed number of excess events is due to a 
fluctuation of the expected background is between 7.8 x 10“® and 1.8 x 10“^, depending on 
different selection criteria applied to the data. 

Some remarks are in order regarding this analysis for the NEV experiments. These 
experiments do not see any evidence for oscillations and hence, they obtain no information 
on a mass-squared difference; only an upper bound on the oscillation amplitude can be 
derived. Therefore, we consider at a fixed value^° of Hence, the Ay^-values 

shown in the table have to be evaluated for 3 d.o.f. (6 *lsnd, dg, in the (2-1-2) scheme and 
for 2 d.o.f. for (3-1-1) because of Eq. (0). Depending on the mass scheme we £x Am^g^D 
at the best fit values given in Tab. ||; for (3-1-0) we use the best £t value of Am^s^D ^^e 
(3-1-1) scheme. Although the NEV experiments are in agreement with no oscillations, the 
value Ay^gy = 2.3 for (3-1-0) shows that a small contribution of can improve the £t 

slightly. 

Table Q conhrms the results of Secs. |VH| and |VH1] . A combination of only solar and 
atmospheric data prefers the (3-1-1) schemes. Therefore, these data are in perfect agreement 
with the global £t in (3-1-1), but £t worse in (2-1-2). On the other hand, atmospheric and 
SBL data prefer the (2-1-2) scheme; the (3-1-1) best fit point is somewhat in disagreement 
with NEV data. This conflict between different data sets does not show up in the y^-values 
given in Eq. (^3[), since the coupling between the data sets is rather weak. As discussed 
in Sec. ^ only the parameter rjs is common to solar and atmospheric oscillations, only the 
parameter links atmospheric and SBL data, while there is no direct coupling of solar and 
SBL data. All the other 7 parameters can be adjusted to give a good fit of the corresponding 
data set. The remaining conflict between the data sets is completely washed out by the large 
number of data points, which are fitted perfectly. 


XI. CONCLUSIONS 

We have presented a unified global analysis of current neutrino oscillation data within 
the framework of four-neutrino mass schemes, paying attention to the inequivalent classes 
of (3-1-1) and (2-1-2) models. We have included all data from solar and atmospheric neutrino 
experiments, as well as information from short-baseline experiments including LSND and 
the null-result oscillation experiments. We have mapped the leptonic mixing matrix into 
a set of parameters in such a way that they have a well-defined physical meaning in each 
data set, independently of the mass scheme ((3-1-1) or (2-1-2)) considered. For example, one 
of these parameters is t]s, the fraction of sterile neutrinos in solar oscillations. Similarly 624 , 

Note that the original analyses of the Bugey and CDHS |Q collaborations were performed in this 
way. 
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describes the fraction of in atmospheric oscillations and sin^ 26 *lsnd is characterizing the 


SBL v;. 


(-) 


Vf. amplitude. The fact that it is possible to describe the results of any of the 
given set of experiments in terms of physical quantities independent of the mass scheme 
implies that none of the considered data sets (solar, atmospheric, SBL appearance or SBL 
disappearance) can be used on its own to discriminate between different mass spectra. This 
follows from the approximation ~ 1, which is motivated by the bounds from reactor 
neutrino experiments, and from the strong hierarchy among the mass-squared differences 
indicated by the data. We have shown how the differences between the mass schemes 
manifest themselves only when two or more data sets are combined. 

We have found that combining only solar and atmospheric neutrino data, the (3-1-1) type 
schemes are preferred, whereas atmospheric data in combination with short-baseline data 
prefers (2-1-2) models. By combining all data in a global analysis the (3-1-1) mass scheme 
gives a slightly better £t than the (2-1-2) case, though all four-neutrino schemes are presently 
acceptable. The LSND result disfavors the three-active neutrino scenario with only Am^QL 
and at 99.9% C.L. with respect to the four-neutrino best £t model. We have also 

performed a detailed analysis of the goodness of fit in order to identify which sub-set of the 
data is in disagreement with the best £t solution in a given mass scheme. 

We have found that, in isolation, the LSND data play a crucial role in suggesting the need 
for a four-neutrino scenario, at odds with every other piece of information. The upcoming 
MiniBooNE experiment will test the very important result of LSND in the near future. 
However, we have shown in this work that the existing data cannot decide between (2-1-2) 
and (3-1-1) mass schemes in a statistically signihcant way. Most likely this problem will 
remain also if MiniBooNE would conhrm the LSND result, and to resolve the ambiguity 
more experimental information will be needed. 

Such information could be provided by experiments with a high sensitivity to the sterile 
component in solar and/or atmospheric neutrino oscillations. One possibility to improve this 
sensitivity for atmospheric neutrinos could be the consideration of neutral current events in 
atmospheric neutrino experiments [Q. Concerning solar neutrinos, we note that the dif¬ 
ferent oscillation solutions show very distinct behavior with respect to a sterile component 
(see Fig. H); hence, the identification of the true solution is important. Moreover, improved 
measurements of neutral current event rates may increase the sensitivity to sterile oscilla¬ 
tions. However, as shown in Ref. |^, the information obtainable from the neutral current 
measurements currently performed at SNO will be limited because of the relatively large 
uncertainty in the flux of solar ®B neutrinos. On the other hand, more data on and/or 

SBL disappearance probabilities could help to solve the (3-1-1) versus (2-1-2) puzzle. Es¬ 
pecially the existing bounds on disappearance are rather weak. This will be improved 
by the MiniBooNE experiment, which - beside testing the Vf. appearance channel - 

will also provide a new measurement for the survival probability ||59|| . 

In view of the ambiguities implied by present data we are looking forward to the results of 
future neutrino oscillation experiments, which may unravel the secret behind the structure 






of the leptonic weak interaction. 
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